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PREFACE 

A  surprising  variety  of  seemingly  unrelated  problems  can  oe 
described  by  the  analog  of  flow  in  a  network.  Analysts  were  presented 
with  both  a  very  general  algorithm  arid  a  very  efficient  code  for  finding 
solutions  to  a  variety  of  network  flow  problems,  with  the  distribution 
of  R.  Clasen's  computer  code  for  this  algorithm  (SHARE,  RS  OKFI) .  and 
the  publication  by  D.  R.  Fulkerson,  "An  Out-of-Kilter  Method  for  Minimal 
Cost  Flow  Problems,1’  J.  Soc.  Indust.  Appl.  Math.,  Voi.  9,  1961,  pp. 

18-27.  In  discussing  the  use  of  two  separate  models  with  various  Air 
Force  organizations,  one  author  has  pointed  out  the  computational  advan¬ 
tages  of  network  flew  models  over  more  general  mathematical  programming 
models,  and  has  become  aware  that  many  operations  analysts  want  to  better 
und  -stand  the  Out-of-Kilter  algorithm;  see  E.  P.  Durbin,  An  Interdiction 
Model  of  Highway  Transportation,  The  RAND  Corporation,  RM-4945-PR,  May- 
1966;  and  E.  P  Durbin  and  Olivia  Wright,  A  Model  for  Estimating  Military 
Personnel  Rotation  Base  Requirements,  The  RAnD  Corporation,  RM-5398-PR, 
October  1967.  This  present  Memorandum  describes  the  operation  and  capa¬ 
bility  of  the  Out-of-Kilter  algorithm  under  the  assumption  that  the 
reader  is  conversant  with  basic  ’inear  programming.  It  should  be  useful 
to  analysts  and  planners  interested  in  a  versatile  modeling  concept  and 
computational  tool. 

D.  M.  Kroenke  is  a  cadet  at  the  United  States  Air  Force  Academy, 
and  coauthored  this  Memorandum  while  on  temporary  assignment  to  The  RAND 
Corporation. 
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SUMMARY 


The  analog  of  steady  state  flow  in  a  network  of  nodes  and  arcs 
may  describe  a  variety  of  processes.  Familiar  examples  are  transpor¬ 
tation  systems  and  personnel  assignment  actions.  Arcs  generally  have 
cost  and  capacity  parameters,  and  a  recurring  problem  is  that  of 
determining  a  minimum  cost  route  between  two  points  in  a  capacitated 
network.  I f  .a  process  can  be  modeled  as  a  network,  and  the  criterion 
for  evaluating  performance  of  the  process  can  be  related  to  the  variables 
corresponding  to  flows  in  the  network,  then  determining  a  minimum  cost 
flow  is  ■equivalent  to  determining  an  optimal  set  of  variables  for  the 
process , 

Efficient  and  general  methods  of  solving  the  minimum  cost  flow 
problem  are  therefore  useful  and  important.  Fu..  .rson's  Out-of-Kilter 
Algorithm  is  an  extremely  efficient  and  general  method  for  solving 
such  problems.  The  algorithm  operates  by  defining  conditions  which 
must  be  satisfied  by  an  optimal  "circulation"  in  a  network  --  roughly, 
a  flow  which  satisfies  capacity  restrictions  on  all  arcs  and  also 
satisfies  stated  conservation  of  flow  conditions  at  all  nodes.  When 
such  an  optimal  circulation  is  determined,  all  arcs  are  "in-kilter." 

At  some  point  in  the  operation  of  the  algorithm,  if  such  a  circulation 
does  not  yet  exist,  some  arcs  are  "out-of-kilter"  --  hence  the  name  of 
the  algorithm.  The  algorithm  arbitrarily  selects  an  out-of-kilter  arc, 
and  tries  to  rearrange  flows  to  bring  that  arc  into  kilter  while  not 
forcing  any  other  arc  farther  out-of-kilter.  If  the  out-of-kilter  arc 
can  be  brought  into  kilter,  '‘he  algorithm  selects  another  out -of-kilter 
arc  and  repeats  fl~i!  process.  Since  there  are  only  a  finite  number  of 
arcs,  repetition  of  this  procedure  eventually  results  in  an  optimal 
solution.  If  any  arc  cannot  be  brought  into  kilter,  the  problem  cannot 
be  solved. 

The  Out-of-Kilter  algorithm  also  solves  the  special  network  problems 
of  finding  maximum  flow  between  two  nodes  in  a  costless,  capacitated 
network,  and  finding  the  shortest  route  between  two  points  in  l  network. 
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I .  NETWORKS 

Steady  state  flows  in  networks  may  represent  many  physical  and 
nonphysical  systems.  The  complex  of  freeways  shown  in  Fig.  1,  linking 
New  York,  Chicago,  St.  Louis,  Denver,  Los  Angeles,  Houston,  Seattle, 
and  Washington,  D.  C.,is  a  network.  The  vehicles  moving  over  that 
network  may  be  considered  homogeneous  units  of  flow.  Alternatively, 
this  network  might  describe  a  petroleum  distribution  system.  Flows  in 
networks  can  also  represent  communications  between  people  in  an  organi¬ 
zation,  inventory  and  production  smoothing  processes  over  time  and  as¬ 
signment  of  personnel  to  jobs. 


Fig.  1  --  A  Network  Linking  Cities 
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The  comnon  elements  of  these  situations,  and  hence  the  abstract 
definition  of  a  network,  are  a  collection  of  points  called  "nodes," 

(the  cities  in  Fig.  1)  and  a  collection  of  "arcs"  wh*ch  connect  these 
nodes  (the  highways  or  pipelines  in  Fig.  1).  We  denote  the  nodes  by 
single  lower  case  letters;  for  example,  node  i,  and  identify  the  arcs 
by  naming  the  nodes  they  connect,  arc  (i,  j) .  Some  homogeneous  com¬ 
modity  (vehicles,  petroleum)  can  flow  over  the  arcs,  and  we  denote  by 
the  amount  of  commodity  flowing  on  arc  (1,  j)  from  node  i  to  node 

).  If  x  ,  <  0  then  the  commodity  flows  from  j  to  i. 

J  lj 

In  most  network  problems,  arcs  have  cost  and  capacity  character¬ 
istics.  Generally,  some  cost  is  incurred  to  move  a  unit  from  node  i 
to  node  j,  and  we  denote  these  unit  movement  costs  by  c^.  This  may 
be  dollars  per  unit  pumped  in  the  petroleum  distribution  network.  We 
also  frequently  find  that  flow  is  limited  by  upper  bounds  or  capac¬ 
ities  on  the  arcs.  For  example,  only  a  limited  number  of  vehicles  per 
hour  can  move  through  the  Lincoln  Tunnel,  and  a  limited  number  of 
barrels  per  day  can  move  througn  the  pip'-Hne  from  Houston  to  St,  Louis, 
We  denote  these  maximum  arc  capacities  u  .  There  may  also  be  a  re- 

*  J 

quirement  for  a  minimum  amount  of  flow  al^ng  any  arc.  We  denote  this 
by  I  .  Imposing  thi6  condition  allows  us  to  construct  networks  with 
controlled  flows  that  describe  particular  problems,  or  that  may  describe 
actual  minimum  demand  levels  at  points  in  a  commodity  flow  network. 

To  summarize,  a  network  is  characterized  by  nodes,  i;  arcs  between 

the  nodes  (i,  J) ;  flow  across  the  arcs,  x^.,  unit  costs  of  flow  across 

the  arcs,  c^;  upper  bounds  on  flow  across  the  arcs  u  ;  and  lower 

bounds  on  flow  acrois  the  arcs,  t  These  characteristics  can  com* 

t  j 

pletelv  characterize  steady  state  flow  in  a  network. 

In  a  problem  with  no  costs,  we  allow  c  -  0.  In  a  problem  uith 
no  lower  bounds  but  flow  only  In  one  direction,  we  allow  f  ■  0.  If 
there  are  no  upper  bounds,  we  allow  u, .  ■  +  ®.  We  further  sasume  that 
all  costs,  flows,  and  bounds  are  integers.  This  is  not  an  overly  strong 
assumption  since  it  all  numbers  are  rational,  clearing  fractions  will 
yield  integers.  For  computat iona 1  purposes  restriction  to  rationals 


suftices.  The  assumption  of  Integra  1 -valued  parameters  is  used  to 
demonstrate  convergence  ot  the  Out-of-Kilter  algorithm. 


-3- 


A  general  network  problem  is  finding  a  minimum  cost  circulat ion 
in  a  network  with  arc  capacities.  In  Sec,  II  we  explain  Fulkerson's 
"Out  -o  f -Ki  1  te  r"  algorithm  developed  to  solve  tins  problem.  The  problem 
requires  that  we  find  flows,  x  ,  that  minimize  total  cost 

( 1 )  j  c  .  .  x .  .  ,  for  all  i  and  j 


while  at  the  same  time  satisfy  the  constraining  conditions 


l  <  x  <  u.  .  for  all  i  and  ]  , 
t  J  i  J  i  J 


and  show  that  in  a  c  i  rcu  la  t  ion ,  wiiat  goes  into  a  node  must  come  out  of 
the  node.  This  is  representtd  by 


0  f  o  r  a  i  1  i , 


While  there  an  several  ways  of  solving  this  problem,  and  if  mav 
be  conceptually  viewed  as  a  linear  programming  problem,  the  Out -of -Kilter 
algorithm  is  both  the  most  general  of  the  specialized  algorithms  and  is 
far  more  efficient  here  than  a  standard  linear  programming  algorithm 
would  he.  We  disc,  ss  some  capabilities  ot  the  Out  -of-Ki  Iter  algorithm 
in  Sec,  III,  and  note  he  re  only  two  reasons  tor  emphasizing  its  ufililv; 

(  1  all  problems  solvable  by  the  more  special  i  zed  network  algorithms 
(maximum  flow  through  a  capacitated  network ,  shortest  route  through  a 
network,  etc.)  are  solvable  by  the  Out -of -Ki Iter  algorithm,  and  an 
understanding  of  it  allows  const  rue  t  ion  ot  the  snore  ef:  icient  specialized 
algorithms  ;  ( 2)  ma  ny  problems  at  e  simple  too  large  to  Ke  hand  led  go 
putat ionally  by  linear  programming  algorithms  on  current  computers . 
Problem  formulation  as  a  networ1  permits  el  tic  tent  cowpt.tal  iorn  l  analysis 
hv  ? he  Out -ol-Ki Iter  algorithm. 
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IT.  THE  OUT-OF-KILTER  ALGORITHM 


INTRODUCTION 

A  set  of  flows,  x  ,,  that  satisfies  flow  constraints  (2)  and 
i  j 

conservation  of  flow  equations  (3)  is  called  feasible.  In  attempting 

to  find  a  minimal  cost  feasible  circulation,  the  Out-ot-Kilter  algorithm 

(OKA)  operates  with  both  arc  costs,  c  ,  and  "r  prices,"  ^ }  which 

arise  from  duality  coca idera t ion6 ,  We  assume  the  reader  has  no  knowledge 

of  duality  theory,  and  attempt  to  make  the  existence  of  the  node  prices 

* 

plausible  and  their  nature  clear  by  means  of  an  economic  explanation. 

We  have  previously  described  the  minimal  cost  circulation  problem  from 
the  point  of  view  of  a  producer  or  distributor.  Consider  the  petroleum 
distribution  network  of  Fig.  1.  Assume  some  regulatory  agency,  say  the 
"Federal  Petroleum  Distribution  Commission,"  ha6  set  this  distributor 
a  variety  of  minimum  and  maximum  flow  levels  on  various  route  segments. 
The  minimum  flow  levels  may  reflect  consumer  d  mands  at  points  in 
the  network,  whereas  the  maximum  flow  levels  may  reflect  physical  rou^e 
capacities.  The  distributor  knows  the  unit  transit  charges  on  each 
route  segment  and  must  decide  the  route  structure  --  arcs  and  the  amounts 
moving  on  each  arc  --  to  meet  the  induced  flow  requi  ements  or  all  arcs 
at  minimum  transportation  cos. 

The  minimum  flow  levels  on  various  route  segments  have  presumably 
been  established  at  the  insistence  of  onsumers.  Therefore  let  o 
denote  the  price  a  consumer  at  node  i  must  pay  for  a  unit  of  petroleum. 
The  will  be  related  to  the  amount  demanded  at  some  nodes,  since  as 
the  amount  demanded  increases,  the  overall  difficulty  in  supplying  it 
will  increase. 

The  Commission  recognizee  that  the  nrice  of  petroleum  to  the  con¬ 
sumer  in  ea-h  city  of  the  network  should  be  related  to  the  distributor's 
overall  cost  in  establishing  the  entire  route  structure,  and  not  re¬ 
lated  to  just  the  transportation  costs  connecting  the  source  and  node  i. 


For  a  further  discussion  of  duality  see  L.  R.  Fv>rd  and  D.  R, 
Fulkerson,  Flows  in  Ne tworks ,  Princeton  University  Press,  Princeton, 

New  Jersey,  1962,  pp.  26-30  and  pp.  113-113.  Also  see  W,  J.  Baumol, 
Economic  Theory  and  Operations  Analysis,  Prentice-Hall,  New  York,  i 965, 

Chap.  6. 
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Neicher  the  distributor  nor  the  Commission  are  initially  sure  what 
these  node  prices,  "  ,  ought  to  be,  but  the  Commission  --  hopefully 
concerned  with  consumer  welfare  as  well  as  producer  profit  --  Intends 
to  see  that  consumer  costs  are  no  greater  than  necessary  to  induce  the 
distributor  to  establish  the  route. 

The  OKA  may  be  pictured  as  a  sequence  of  systematic  decisions  the 
distributor  makes  with  the  Distribution  Commission  watching  carefully. 

At  each  step  the  distributor  asks,  "'..'hich  route  shall  be  used  to  mini- 
mire  total  system  costs,  taking  account  not  only  of  transportation  costs, 
but  also  of  the  commodity  prices  at  the  various  cities?"  If  there  is 
no  profitable  route  at  some  step,  the  distributor  tells  the  Commission, 

’’You  must  now  allow  increased  prices  at  some  cities  in  order  to  make  it 
financially  possible  for  me  to  continue  constructing  a  route."  The 

3 

Commission  answers,  "Okay,  but  let's  find  that  route  which  requires  the 
minimum  price  increase." 

Based  on  this  fanciful  exchange,  we  associate  with  each  node,  i, 
a  variable,  n  that  can  be  considered  the  price  of  a  unit,  of  the  flow 
commodity  at  the  node.  The  Distribution  Commission  defines  a  net  arc 
cost,  c  ,  as 


(A) 


‘ij 


+  TT 


n.. 

J 


The  new  cost,  c  ,  represents  the  total  cost  to  the  system  -- 
consumer  and  distributor  --  of  transporting  one  unit  of  flow  from  node 
i  to  node  j.  This  definition  compares  the  cost  of  retaining  a  unit  at 
node  i  with  the  cost  of  moving  it  to  node  j.  In  moving  a  unit  of  flow 
from  i  to  j ,  the  commodity  price  at  i. ,  T'^  >  is  foregone,  and  an  actual 

transportation  cost,  c  ,  incurred.  If  the  sum  of  these  costs  is 

greater  than  the  commodity  price  at  j ,  ,  then  it  does  not  pay  to  ship 

a  unit  from  1  to  J.  The  c  will  be  positive.  On  the  other  hand,  if 
a  unit  at  j  costs  more  than  at  i  plus  the  transportation  cost,  c  will 
be  negative  --  the  system  benefits  from  the  move  --  and  shipment  from 
i  to  j  is  profitable.  If  the  value  at  j ,  n ( ,  is  balanced  exactly  by 

J 

the  value  at  i  plus  the  transportation  cost,  (tt^  -f  then  c  *= 

and  we  are  indifferent  to  an  additional  unit  flowing  from  i  to  j. 


0 
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Limitatlons  on  permissible  flow  levels  (2)  together  with  possible 
levels  of  total  system  cost  (A)  yield  the  following  conditions  that 
will  be  satisfied  by  an  optimal  solution  to  the  minimal  cost  circulation 
problem  (1)  -  (3). 


(5) 

If 

c  < 

ij 

0 

then 

X’J 

"  uir 

(6) 

If 

Cij  ” 

0 

then 

9 

“ij 

-xij-uir 

(7) 

If 

V 

0 

then 

Xij 

Equation  (5)  states  that  when  net  arc  cost  is  negative  (when  it  is 
profitable  to  send  the  commodity  from  1  to  j) ,  flow  on  the  arc  ought  to 
be  as  large  as  possible.  Equation  (6)  states  that  when  net  arc  cost  is 
xero,  we  arc  indifferent  to  flow  level  so  long  as  it  meets  constraints. 
Equation  (7)  states  that  when  net  arc  cost  is  positive  (a  loss  is  in¬ 
curred  „y  sending  the  commodity  from  i  to  j) ,  flow  on  the  arc  ought  to 
be  at  the  minimum  level  possible.  The  OKA  is  designed  to  construct  a 
flow  meeting  these  conditions. 

Any  arc  that  meets  the  optimality  conditions  (5),  (6) ,  (7)  is 
defined  as  "in-kilter,"  Arcs  that  do  not  satisfy  these  conditions  are 
"out-of-kilter." 

Out -of -ki Iter  arcs  can  be  grouped  into  two  categories: 

(a)  Those  that  are  feasible  but  not  optimal.  They  have  flow  which 
satisfies  (2),  but  prices  and  flow  do  not  satisfy  (5),  (6),  (7). 

(b)  Those  that  are  infeasible.  Flow  is  below  or  above  the  upper 
bound.  Therefore  (2)  is  not  satisfied. 

Arcs  that  are  feasible  but  not  optimal  must  fit  one  of  the  follow- 

ie 

ing  states  or  conditions. 


It  is  not  necessary  to  consider  as  a  special  case  *  <  0  and 
xij  <  ^ij  or  Cjj  >  0  and  x^,  >  uij-  These  cases,  though  infeasible, 
are  included  in  conditions  I  and  II,  If  the  algorithm  terminates  with 
all  arcs  in  kilter,  any  originally  infeasible  arc  in  state  I  or  II  will 
not  only  be  feasible  but  optimal. 
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Conditiun  I; 
Condition  II: 


Condition  III: 


Condition  IV: 


Condition  V: 


Condition  VI: 


<  0 

and 

< 

ij 

ij 

ij 

c,  , 

O 

A 

and 

> 

£ 

ij 

ij 

ij 

one  o 

f  these  conditions: 

c ,  , 

>  0 

anc 

X  ,  . 

< 

£  .. 

ij 

ij 

ij 

c ,  . 

1 

o 

and 

x  r  , 

< 

l ,  , . 

IJ 

ij 

ij 

c ,  , 

*•  0 

and 

x  . 

> 

u ,  ,  . 

ij 

ij 

ij 

c  , 

<  0 

and 

x ,  , 

> 

u ,  , . 

ij 

ij 

ij 

A  "kilter-number,"  as  defined  below,  is  associated  with  each  arc. 


Arc 

Condition  Kilter  Number 

I  . •  “i)] 

II  . ctjLXij  "  ‘lj’* 

Ill,  IV .  [liJ  -  x  ] 

V«  VI .  [xij  '  ui^ 


Notice  that  in  all  cases,  the  kiltev  number  is  positive.  For  the 
feasible  states,  I  and  II,  the  kilter  number  is  a  measure  of  non¬ 
optimality,  The  kilter  number  for  states  III  -  VI  indicates  the  degree 
to  which  an  arc  is  infeasible.  In-kilter  arcs  meet  conditions  (5), 

(6),  or  (7),  and  therefore  have  a  kilter  number  of  zero. 

The  OKA  operates  by  arbitrarily  selecting  an  out-of-kilter  arc  and 
rearranging  flows  in  an  attempt  to  r.duce  the  killer  number  of  that  arc 
to  zero.  During  this  process  the  kilter  numbers  of  other  arcs  do  not 
increase,  and  may  in  fact  decrease.  The  algorithm  terminates  with  an 
optimal  solo' ion  when  kilter  numbers  of  all  arcs  are  zero.  Optimality 
of  such  a  solution  is  verified  by  recalling  the  definition  of  c^,  and 
noting  that  a  kilter  number  of  zero  Implies  Eqs.  (5),  (6),  and  (7)  hold. 
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It  is  slightly  more  complicated  to  show  that  the  algorithm  termi¬ 
nates.  The  algorithm  selects  one  out-of-kilter  arc  at  a  time  and 
attempts  to  bring  that  arc  into  kilter  by  rearranging  flows  without 
forcing  another  arc  out  of  kilter.  If  this  process  can  be  completed, 
the  algorithm  seeks  another  out-of-kilter  arc.  If  it;  cannot  be  com¬ 
pleted,  the  problem  is  infeasible.  If  it  takes  only  a  finite  number 
of  steps  to  put  an  arc  into  kilter  or  to  determine  infeasibility,  and 
there  are  only  a  finite  number  of  arcs  in  the  network,  the  algorithm 
must  terminate.  We  will  later  indicate  that  these  conditions  are 
satisfied, 

OPERATION  OF  THE  ALGORITHM 

Once  the  minimal  cost  circulation  problem  (2)  has  been  formulated, 
the  OKA  can  be  started  with  any  set  of  node  prices,  ,  and  any  circu¬ 
lation  which  satisfies  the  conservation  of  flow  equations  (3).  This 
circulation  and  set  of  node  prices  can  initially  be  zero.  Arbitrarily 
select  an  out-of-kilter  arc  joining  two  nodes,  say  s  and  t.  The  fact 
that  the  arc  (s ,  t)  is  out  of  kilter  indicates  that  it  is  either  profit¬ 
able  or  necessary  (or  possibly  both)  to  ship  an  additional  unit  from 
t  to  s  or  from  s  to  t .  In  either  event,  flow  change  is  always  desired 
for  an  out-of-kilter  arc.  In  order  to  change  flow  on  the  arc  and  yet  keep 
flow  in  the  network  balanced  (3),  another  path  through  the  network  from 
t  to  s  must  be  found  along  which  flow  values,  x  ,  can  be  changed.  In 
conetructing  this  path,  tv.  pieces  of  information  (a  "label")  must  be 
retained  at  every  node.  The  first  component  of  the  label  at  a  given 
node  )  indicates  the  previous  node  in  the  path,  and  whetntr  present  flow 
moves  from  i  to  j,  denoted  (i+),  or  from  j  to  i,  (i-).  The  second  com¬ 
ponent  of  the  label  is  the  amount  b,  which  flow  on  arc  (i,  j)  is  to  be 
changed,  *{)).  A  complete  lab*l  at  node  j  is 

[ i± ,  c(j)l. 


The  search  proceeds  from  node  t  through  the  network  seeking  a  path  back 
to  s.  "labels"  are  assigned  to  various  nodes  so  that  when  a  path  is 
found  connecting  t  and  s,  the  labels  indicate  the  direction  and  magni- 
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*-ude  of  necessary  flow  change  along  the  path.  Of  course  not  every 
>de  will  receive  «  label  in  this  process.  Circulation  from  Houston  to 
St.  Louis  cstn  be  completed  by  arcs  from  St.  Louis  to  Washington  and 
Washington  to  Houston.  Seattle  and  Los  Angeles,  for  example,  would  not 
receive  labels. 

If  a  path  connect itig  t  end  s  is  not  found,  the  algorithm  determines 
new  node  prices,  tr  in  sue  a  manner  that  (a)  another  node  will  be 
labeled  in  the  partial  (f ,  e)  path,  (b)  one  le6S  arc  will  be  con¬ 
sidered  for  inclub  ion  in  the  (t,  s)  path,  or  (c)  if  no  arcs  remain  and 
the  path  i6  incomplete,  the  problem  is  deemed  infeasible. 

We  next  explain  some  of  the  labeling  rules  that  are  used  in  at- 

■k 

tempting  to  find  a  path  in  the  network  from  t  to  s .  If  a  path  is 

found  we  term  this  "breakthrough."  Assume  that  an  arbitrary  arc  (s ,  t) 

is  in  state  I,  with  c  <0  and  x  <  u  Since  c  <  0  it  is  profit- 

st  st  8 1  st 

able  to  increase  the  flow,  ,x  .  We  begin  at  t  by  assigning  node  t  the 

S  L 

label  [s+,  *(t)l  where 

(8)  «(t)  "  V  -  xst. 


This  indicates  that  flow  on  this  arc  should  be  increased  by  e(t)  in 
order  to  be  equal  to  : he  upper  bound,  and  hence  be  optimal.  See  Fig.  2. 


's+,  € ( t ) ] 


Fig  2  --  Initial  Flow  on  Out-of-Kilter  Arc  (s,  t) 


We  describe  cases  I  and  IV.  The  rules  are  onlv  slightly  different 
in  the  other  cases,  and  the  method  remains  basically  the  same.  The  im¬ 
portant  special  network  flow  problem  which  case  IV  covers  is  the  "capaci¬ 
tated  transportation  problem."  Section  III  describes  operation  of  the 
algorithm  in  this  problem. 
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The  algorithm  next  seeks  an  unlabeled  node  adjacent  to  t,  say  j. 
Assume  conditions  on  the  arc  (c,  J)  are  either 


(9a) 


c  ,  >  0  and  x  ,  <  ( 
tj  tj  tj 


or 


(9b)  c  £  0  and  x  .  <  u  .. 

M  tj  tj 

In  either  case  we  desire  to  increase  flow  enough  from  t  to  j  to  correct 
infeasibility  or  nonoptimality  on  the  arc.  We  therefore  assign  node  j 
the  label  [t+,  e(j)\  If  condition  (9a)  holds,  we  have 

dA'0  e(.j)  x  min[e(t)  ,  £  -  x,.  ] , 

while  for  condition  (9b)  we  have 


(10b) 


e( j)  *  min[e(t) ,  u  -  x  ' 


The  e(j)  state  the  maximum  amount  by  wnich  we  can  raise  flow  on  the 

arc  (t .  j) •  See  Fig.  3.  We  would  like  to  increase  x  bv  ( £  -  x  ) 

-  tj  tj  tj' 

in  case  9a.  And  we  would  like  to  increase  x  by  (u  -  x  )  in  case 

^  J  ^  J  t  j 

9b.  But  at  the  same  time  the  previous  label  at  t  states  that  we  want 
to  increase  flow  on  (s ,  t)  at  most  b\  the  amount  e  ( t )  .  Thus  we  compare 
tvO  with  the  new  value  computed  by  (10a)  or  (10b)  and  select  e(j)  as 
the  mlnlmui  of  the  two  values  --  the  maximum  amount  by  which  we  can  in¬ 
crease  flow  along  both  arcs  without  violating  the  feasible  bounds  or 
optimality  conditions  on  any  arc.  This  guarantees  that  flow  changes 
will  not  increase  the  kilter  number  ot  anv  arc. 


Fig.  3 


Carrying  Forward  a  Label 
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Assume  that  we  next  encounter  an  arc  connecting  node*  s  and  j.  We 
have  then  completed  a  path  from  e  co  t  to  j  and  on  returning  to  a  "long 
the  path  we  wish  to  increase  flow.  Thus  far  the  flows  encountered  have 
been  in  the  appropriate  direction--s  to  t  and  t  to  j.  Suppose,  however, 
that  flow  on  the  arc  (s,  J)  is  moving  from  s  to  j ,  and  that  arc  (s ,  j) 
is  in  one  of  the  following  conditions: 


(11a) 

(lib) 


,  2  Cl 

and 

X  >  J t 

s  j 

sj 

SJ 

.  <  o 

and 

x  , 

>  u  , 

sj 

Sj 

Sj 

We  then  assign  the  label  [j-,  €(*)]  where 


(12a)  e(s)  *  min[*(J),  x  -  t  if  (11a)  holds 

®J  s.j 


or 


(12b) 


e(s)  -  mini *(j)  ,  x  -  u  ' 


if  (lib)  holds. 


In  either  case  we  want  to  increase  flow  on  the  directed  path  s, 

(s,  t)  ;  t,  (t,  j)  ;  J,  (J,  c)  by  the  amount  e(s).  We  therefore  decrease 
flow,  x  ,  by  the  amount  e(s),  This  is  equivalent  to  increas ing  iiow 
in  the  desired  direction  (see  Fig.  4)  by  e(s).  Recall  that  c(s)  is 
the  >.*aximun>  amount  by  which  flow  or  all  arcs  can  be  increased  without 
increasing  any  kilter  number. 

We  have  now  broken  through  --  found  a  path  from  t  to  s  through  the 
network ,  and  have  determined  e(s)  ,  the  maximum  amount  by  which  flow 


* 

Consider  the  petr.  .pirn  distribution  network  of  Fig.  1.  Assume 
present  flow  is  6  units  from  Houston  to  Denver  and  4  units  from  Houston 
to  Washington,  D.C.  It  we  desire  to  increase  flow  from  Denver  to 
Houston  by  2  units,  it  is  rot  reasonable  to  send  b  units  f ro.m  Houston 
to  Denver  and  at  the  same  time  to  send  2  units  of  the  same  commodity 
from  Denver  to  Houston.  It  would  be  more  profitable  to  reduce  the 
i low  from  Houston  to  Denver  to  4  units.  This  "increases"  the  flow 
from  Denver  to  Houston  by  reducing  the  amount  Houston  must  ship  to 
Denver,  and  is  an  example  of  increasing  forward  flow  by  decreasing 
reverse  flow. 
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FJfe.  4  --  Labeling  an  Arc  with  Reverse  Flow 

along  this  path  can  be  increased  without  violating  optimality  or  feasi¬ 
bility  conditions.  The  algorithm  now  changes  flow  according  to  the 
label  at  each  node.  At  t  we  increase  flow  from  s  to  t  by  e(s).  At 
j  we  increase  flow  from  t  to  j  by  e(s).  And  at  s  we  dec reabe  flow 
from  s  to  j  by  e(s).  Flow  from  s  to  t  has  now  increased  by  £,.)•  The 
amount  of  flow  augmentation,  e(s)  ,  was  chosen  bv  considering  for  each 
arc  the  amount  by  which  flow  could  change  without  violating  feasibility 
or  optimality  conditions.  Therefore,  on  at  least  one  arc,  the  one 
yielding  the  minimum  of  the  e(i),  flow  has  actually  been  set  equal  to 
an  upper  or  lower  bound.  Reviewing  the  definition  of  the  kilter 
numbers  for  conditions  I  -VI,  and  the  flow  augmenting  rules,  note 
th„:  flow  augmentation  makes  flow  on  at  least  one  arc  closer  to  a  bound , 
and  therefore  at  least  one  kilter  number  has  decreased. 

We  summarize  operation  of  the  algorithm  to  this  point.  After 
starting  with  any  circulation  and  node  prices,  we  arbitrarily  select 
an  on  t -o  f -k  i  1 1  er  arc  y  -•  >  O  and  label  it  t.  We  next  examine  all  nodes 
ad jacenC  to  t  ,  determine  if  1  low  on  arcs  to  these  nodes  could  he  ap¬ 
propriately  modified  to  help  bring  (s,  t)  into  kilter,  and  label  the 
nodes  using  rules  9  through  12.  For  each  node  so  labeled  we  attempt 
to  repeat  this  labeling  process  until  a  path  returning  to  s  it  found. 
Flow  along  that  path  is  then  augmented  bv  the  maximum  feasible  amount  , 
€(•>)  ,  which  is  the  minimum  of  feasible  amounts  along  all  arcs  on  the  pat! 

★ 

The  kilter  number  will  not  necessarily  become  zero  on  any  new  arc. 
If  the  minimum  e(i)  occurs  on  an  arc  original ly  satisfying  Kq.  (h), 
with  original  kilter  number  zero,  the  ki  Lter  'lumber  of  that  arc  will 
remain  zero,  and  kilter  numbers  on  other  arcs  will  decrease,  but  not 
necessarily  to  zero.  Notice  also  that  since  flow  ctianges  by  an  in¬ 
tegral  value  and  since  all  costs  and  prices  are  integers,  kilter 
numbers  decrease  by  integral  amounts. 
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All  labels  are  now  erased,  and  arc  (s ,  t)  is  examined  again.  If 
that  arc  is  in  kilter,  we  seek  another  out -of -ki iter  arc.  If  (s,  t) 
remains  out  of  kilter,  the  search  and  label  process  is  repeated.  The 
same  path  from  t  to  s  will  not  be  chosen  again  because  on  at  least  one 
arc  of  that  path,  flow  is  at  an  upper  or  lower  bound.  If  all  arcs  are 
in  kilter,  the  optimality  and  feasibility  conditions  (5),  (6),  and  (7) 
are  satisfied,  and  the  algorithm  has  terminated  with  an  optimal  solution 
to  (1)  ,  (2)  ,  and  (3) . 

Labeling  rules  for  out -o f -ki 1 te r  arcs  (i,  j)  in  states  other  than 
I  or  IV  are  similar.  When  flow  on  arc  (f,  j)  should  be  decreased  (states 
II,  V,  VI),  the  labeling  process  starts  at  i,  attempting  to  move  along 
a  path  through  the  network  to  j.  An  out -of -ki Iter  arc  in  state  III  is 
treated  as  one  state  I  or  IV,  except  that  e(j)  *  i .  .  -  x, . ,  the  imount 

iJ 

of  flow  augmentation  needed  to  remove  the  infeasibility. 

NON  -  BREAKTHROUGH 


Given  an  out-of-kilter  arc  (s,  t)  on  which  we  want  to  increase  flow, 

out  ability  to  find  a  path  through  the  network  from  t  to  s  has  depended 

on  the  availability  of  arcs  (i,  j)  in  one  of  two  conditions.  On  an  arc 

(i,  i)  either  we  both  desire  to  increase  flow  x  and  are  able  to,  or 

1 J 

we  both  want  to  decrease  f’ow  x..  and  are  able  to.  If,  in  attempting 

to  find  a  path  for  a  given  ou t -of  - ki 1 1 er  arc  (s ,  t),  we  have  labeled  all 

nodes  we  can  bv  rules  (9a),  (9b),  (11a),  (.ilb),  and  have  not  returned 

to  node  s,  we  have  a  "non-breakthrough . 

There  are  two  possible  difficulties.  Either  (1)  the  net  arc  costs 

c  provide  insufficient  incentive  for  a  distributor  to  move  the  com  • 

modify  through  the  network,  or  (2)  no  augmentation  is  possible  because 

flows  x  are  equal  to  upper  bounds  and  flows  x  are  equal  to  lower 
i  ]  5 1 

bounds.  If  case  1  occurs,  we  will  s  ys  t  etna  t  i  ca  1 1  y  change  prices  ”  in 
a  way  to  be  described,  recompute  c  ,  and  again  t  rv  to  break  through. 

1  j 

If  case  2  occurs  w<  must  declare  the  problem  iji  feasible.  This  in  - 
leas  ibt  lit  v  ,ar.  he  seen  since  we  reached  this  point  b\  attempting  to 
bring  an  out  -o  I-ki  1  ter  arc  (s,  t )  into  kilter.  The  fact  that  the  rc 
( t>  ,  t a  s  a  positive  kilter  number  indicates  that  either  optimality 
conditions  or  feasibility  conditions  are  not  satisfied  on  it.  If  the 
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problem  is  feasible,  and  only  the  optimality  conditions  fail  to  be 
satisfied,  the  repricing  of  the  n  will  eventually  bring  c  to  zero, 

1  8  C 

and  hence  bring  the  original  arc  into  kilter.  The  problem  is  in¬ 
feasible,  however,  if  flow  x  is  outside  the  bounds  (£  ,  u  ),  and  no 
path  through  the  network  exists  to  change  that  flow  and  bring  x 

s  t 

within  bounds. 

To  look  more  carefully  at  this  non-breakthrough  situation,  con¬ 
sider  the  various  arcs  that  can  be  used  for  a  breakthrough  from  t 
back  to  s.  Arcs  connecting  labeled  nodes  to  labeled  nodes  are 
already  in  the  path  from  t  to  our  present  position  They  cannot  be  of 
further  use.  Arcs  which  connect  unlabeled  nodes  to  unlabeled  nodes  are 
of  no  use,  since  we  have  no  means  of  extending  our  path  to  these  arcs. 

Therefore,  we  need  only  consider  arcs  that  connect  labeled  to  unlabeled 

* 

nodes . 


This  iteration  started  with  an  out-of-kilter  arc  (s,  t)  in  condi¬ 
tion  I  or  IV,  and  we  have  beer  attempting  to  move  flow  fo rwa rd  from  t 
through  the  network  to  6,  and  have  stopped  at  a  labeled  node,  i.  If 
arc  (i,  J)  is  to  be  potentially  useful,  its  flow  must  be  less  than 
capacity  (x^  <  u^j  if  x_  "•  01  ,  or  if  flow  is  moving  in  the  direction 
opposite  to  our  path,  that  flow  must  be  greater  than  the  lower  bound 
(  x  .  .  >  /  , , ) .  In  the  former  case  we  can  increase  forward  t  low ,  and  m 

Ji  Ji 

the  latter  case  we  can  decrease  reverse  flow.  In  either  case  a  useful 
arc  Is  one  that  increases  net  forward  tlow.  For  to  rwa  rd  flow  arcs 

(i,  n  with  x  <  u  ,  arc.  for  reverse  flow  arcs  with,  x  •  i  con- 
1  j  1 1  U  I  i 

<J it  ions  { 9a)  .  (dbl,  (  l  la  1  ,  (lib's  cannot  be  satisfied.  Otherwise  ,  he¬ 
roics  10  and  12  the  notes  j  would  have  received  Libels.  In  other  words, 
the  current  cost  structure  is  incorrect. 


Arcs  witi'.  forward  flow  below  capacit 


ret  arc  cost  (c  ,  01. 

i  j 


Otherwise,  bv  ru 


labeled.  Su-'h  arcs  -.must  also  have  x  o 
would  be  labeled.  Arcs  with  reverse  flow 


v  o.  ,  <  u .  1  have 

1 J  i] 

le  i  9b  ,  bot  h  nodes 

f  .  ,  or  bv  t. da  i  bot 
1 J 

above  the  lower  bo 


pos  i  t  ive 
would  he 
h  nodes 
end 


★ 

If  there  art  no  such  arts  and  t  he  lower  noi:nt!  op  is,  1 1  is  posi¬ 
tive,  the  problem  is  infeasible.  For  s  flows  into  an...  t  flows  to 
the  r  network  nodes,  hut  no  network  node  ever  flows  back  to  s. 
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(x  ,  ?  t  )  must  have  negative  net  arc  coat  (c  <  0)  because  of 
j 1  J *■  J 1 

rule  (11a),  and  must  have  x  £  u  because  of  rme  (lib).  Thus,  the 

arc  conditions  that  must  exist  on  arcs  connecting  labeled  nodes,  i, 

to  unlabeled  nodes,  j,  in  case  of  a  non -breakthrough  are 


(13) 

(14) 


Cij  >  °  3nd  *ij  S  Xij  *  Uij’ 


Eji  <  °-  s  xu  *  Uu- 


Ji  )i 


Equations  (13)  and  (141  represent  arcs  that  nay  yet  be  useful  in 
breaking  through  to  th  stare ing  out-of-kilter  arc  (s,  t).  The  arcs 
are  useful  since  it  is  possible  that  flow  can  stil1  be  modified  on  them. 


Recall  that  Eq.  4  defines  c. 


ij 


ij 


c  +  — 

ij  i 


The  potentially  useful  forward  arcs  have  not  been  used  since  c  >0. 

i  j 

This  states  that  the  value  given  up  at  node  i.  **  ,  plus  the  cost  of 
moving  an  additional  unit  of  f  low  to  i,  c  ,  j  •  ’•eater  than  ~  ,  the 

i  j  j 

value  of  having  the  unit  at  j.  It  does  not  pay  the  distributor  to  send 
it.  For  a  reverse  arc  satisfying  (14)  it  does  not  pay  to  reduce  back  flow. 

At  this  point,  however,  path  construction  has  halted,  and  there 
will  be  no  path,  through  the  network  from  t  to  s  unless  it  been  ,  possi¬ 
ble  to  use  a  previously  unused  arc  satisfying  ( 13)  or  (  1-.)  ,  and  to 
label  at  least  one  previously  ur.lahelec  node.  will  accomplish,  this 

by  raising  the  price  ot  the  commodity,  at  ail  un  labeled  nodes,  j. 

This  will  make  it  mur  profitable  to  ship  toward  the  on  labeled  node  on 
forward  arcs  ( increase  x  )  or  to  reduce  flow  away  :  r  the  uni  a  be  led 


'de  on  reverse  arcs  (reduce  x  )  .  This  repricing  act  i* 


cor  re  a  ponds 

i  x 

to  a  grudging  approval  bv  the  hypothetical  Federal  Petroleum  Mst »  ibut i 
‘ntrar.i  «st  on  in:  an  t  ;u  tease  sn  *  ossrod  i  t  v  prices  that  c  ons-am*  •  s  will  ha.e 
warning  !  r  o-  tin  do-mission:  No  mote  •  •  -so.  ni\  cs.sa: %  to 
i  again.'  We  obviously  cn  :iu  get  1  tow  :;vvu;t  again  m  at 
•v!  at,  s  a  t  v  “tying  -.13)  o:  (  1 4)  bv  adding  a  ve :  v  I  a  t  ge 
:n  1  a  be  :  •  node  prices, 
the  minimum  suf i ic  tent  in 


,  n.e  i.omm ’.  s  s  ion  vi  ,  4  not  a  ;  *ow 


u  .:t  labeled  nod*  prices. 


*  ! 
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We  can  get  flow  moving  again  on  at  ..east  one  appropriate  forward 

arc  by  finding  the  minimum  c  on  all  the  forward  arcs,  say  b, ,  and 

^  J  1 

adding  this  to  ",  at  all  unlabeled  nodes.  This  would  result  in  a  new 
system  cost  on  each  forward  ere: 


ij 


cu  +  h 


<‘j  +b.> 


And  on  the  arc  that  yielded  the  mini 


c  .  -  (c.  . 
ij  iJ 


' .)  -  c  .  . 
j  i  J 


By  labeling  rule  (9b)  this  last  arc  can  be  used,  since  at  system  cost 
of  zero  we  are  indifferent  to  changes  in  flow.  Therefore,  node  j  can 
be  labeled.  This  action  has  essentially  cost  the  consumer  and  the  system 
b,  per  unit  of  flow  arriving  at  j,  which  we  have  had  to  add  to  the  price, 
n  ,  at  the  unlabeled  nodes. 

It  might  have  proven  cheaper  to  reduce  back  flow  on  one  of  the  arcs 

(j,  1).  This  could  have  beer,  accompi  ished  by  finding  b^  ,  the  magnitude 

of  the  smallest  negative  cost  (minjc..j),  and  adding  b  to  all  node  prices, 

j  1  l 

”  .  ,  at  unlabeled  nodes.  On  each  relevant  arc  with  reverse  flow  this 
J 

would  ha.ve  resulted  in  a  new  system  cost  given  by 


(151 


cji 


C  .  .  +  *T  , 

’1  J 


-  ("i  +b2). 


On  the  arc  yielding  ,  the  minimum  magnitude, 


(U) 


ji 


c. . . 

Ji 


=  0. 


By  labeling  rule  (11a),  this  last  arc  can  be  used  to  augment  flow  and 
node  j  can  be  labeled.  The  pricing  procedure  has  essentially  cost  the 
consumer  and  the  system  b,_,  per  unit.  If  b^  is  less  than  b^  we  should 
use  the  second  procedure.  In  general,  then,  the  algorithm  finds  b^  and 
b  from  the  arcs  defined  by  (13)  and  (Id),  selects  b  (the  smaller  of 
h,  and  b^) ,  and  increases  node  prices  at  all  unlabeled  nodes  bv  an 
amount  b. 


Repricing  actions  have  allowed  the  distributor  to  add  another  arc 


to  the  (t  ,  s)  path  at  sor,:e  cost  to  the  purchasers  but  at  the  min inir.n 
cost  to  those  purchasers .  Following  the-  repricing  procedure ,  tie 
algorith:  .  cks  the  original  out-of-kilter  arc  (s,  t)  .  If  it  if  still 
out  of  kilter,  the  search  for  a  path  back  to  node  t  is  repeated.  We 
sunir-arize  the  algorithm  to  this  point.  Arbitrarily  select  an  out-of¬ 
kilter  arc  (s,  t),  and  by  means  of  the  labeling  procedure  attempt  to 
find  a  path  through  the  network  from  t  to  s .  If  such  a  path  is  found, 
appropriately  modify  the  flew  along  each  arc  of  the  path  by  the  amount 
determined  in  the  labeling  process.  If  a  path  to  s  is  not  immediately 
found,  consider  all  arcs  that  are  potentially  useful  in  extending  the 
oath  from  t.  Modify  the  node  prices  on  these  arcs  so  that  one  of  these 
is  used  to  extend  the  path  and  continue  the  labeling  procedure.  If  there 
are  no  arcs  connecting  labeled  to  unlabeled  nodes  that  satisfy  the  con¬ 
ditions  (13)  or  (14),  then  no  path  is  possible  and  the  problem  is  in¬ 
feasible.  Once  this  particular  arc  is  brought  into  kilter,  the  algorithm 
continues  until  all  arcs  are  in  kilter  (which  signifies  an  optimal  so¬ 
lution  to  the  original  problem)  or  until  it  is  impossible  to  construct 
a  feasible  circulation. 

TERMINATION  AND  OPTIMALITY 

The  OKA  concentrates  on  one  out-of-kilter  arc  at  a  time.  Each  time 
a  breakthrough  to  the  out -o f -kilter  arc  occurs,  the  labeling  rules  in¬ 
sure  improvement  in  kilter  numbers  on  every  arc  in  .he  path,  including 
the  out-of-kilter  arc.  We  have  mentioned  the  assumption  that  all  values 
are  integral.  Kilter  numbers  therefore  improve  by  integral  amounts  at 
each  breakthrough,  Th"s  ,  only  a  finite  number  of  breakthroughs  are  pos¬ 
sible  before  the  arc  is  in  kilter.  There  are  only  a  finite  number  of 
arcs.  If  the  algorithm  is  successful  in  bringing  every  arc  into  kilter, 
the  solution  is  optimal  by  (5),  (6),  and  (1). 

The  remaining  situations  are  the  non-breakthroughs.  When  a  non¬ 
breakthrough  occurs  and  potentially  useful  arcs  exist,  price  changes  are 
made  and  at  least  one  of  the  potentially  useful  arcs  is  used  to  extend 
r he  path.  Since  there  are  only  a  finite  number  of  arcs,  we  must 
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eventually  either  break  through  or  discover  there  are  no  potentially 
useful  arcs  left.  In  the  former  case,  we  decrease  the  kilter  number 
cf  the  original  out-of-kilter  arc;  in  the  latter,  we  declare  the 
problem  infeasible. 


Ill,  USES  AND  SPECIAL  CASES 


While  networks  can  be  used  to  model  a  variety  of  actual  problems, 
ingenuity  is  often  called  for  in  formulating  the  network  to  describe 
the  problem.  If  the  network  can  be  properly  formulated,  however,  it  is 
far  mere  efficient  to  solve  a  minimal  cost  circulation  problem  than  the 
'■  juivalent  linear  programming  problem.  Furthermore,  we  frequently  need 
to  examine  the  behavior  of  A  solution  as  the  parameters  vary.  As 

k 

Fulkerson  points  out,  the  Out'Of-Kilter  algorithm  is  designed  to  start 
with  any  circulation  and  any  set  of  node  prices.  Therefore,  a  previously 
derived  optimal  solution  can  be  used  to  begin  a  new  problem  with  re¬ 
sultant  savings  in  computation  time. 

We  have  described  the  OKA  for  cases  I  aid  IV.  The  important  special 

network  flow  problem  case  IV  covers  is  the  "capacitated  transportation 

problem,"  or  the  chipnent  of  a  fixed  level  of  flow,  v  (possibly  maximum 

feasible  flow) ,  through  a  network  from  node  s  to  node  t  at  minimum  cost. 

In  this  special  case  we  can  start  with  all  flows  equal  to  zero,  and  all 

node  prices  equal  to  zero.  On  a  "return"  arc  (t,  s)  we  set  the  lower 

bound,  f,  equal  to  the  upper  bound,  u,  equal  to  the  desired  flow  value, 

v.  Then  the  arc  (t ,  s)  is  out  of  kilter  since  the  initial  flow,  x  -  0, 

is  less  than  the  required  flow,  l  -  v.  The  cost,  c  ,  mav  e.  >itrariiy 

ts  ts  3 

be  set  to  -®,  causing  arc  (t ,  a)  to  be  initially  out  of  kilter  in  state 
I,  or  c  may  be  set  to  0,  causing  arc  (t,  s)  to  be  initially  out  of 
kilter  in  state  IV. 

We  mention  two  other  important  special  cases  of  the  general  minimal 
cost  circulation  problem.  The  first  is  determining  maximum  flow  in  a 
capacitated  network,  and  the  second  is  finding  the  shortest  route  through 
a  network  iri  which  costs  on  arcs  are  times  or  distances.  While  there 
are  specialized  algorithms  for  each  case,  the  Ou t -of -Ki 1  ter  algorithm 
handles  either,  and  in  the  process  indicates  how  to  construct  a  mote 
specialized  algorithm. 

In  order  to  use  the  OKA  to  find  the  maximum  level  of  flow  in  a 
network  from  node  s  to  node  t,  we  take  the  following  actions: 
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(1)  add  an  arc  connecting  t  to  s. 


(2)  set  l. 


(3)  set  c 


ij 

ij 


0  for  all  i  and  \. 

0  for  all  i  and  i,  except  c 


ts 


(4)  set  c  •  -M  where  M  is  very  large, 

L  S 


(3)  set  x. 


ij 


0  and  t,  *  0  for  all  i  and  j. 

l  J 


(6)  set  u  *  M  where  M  is  very  large, 
ts 

On  all  arcs  except  (t,  e)  ,  costs  are  zero,  node  prices  are  zero,  and 
flow  is  feasible  and  zero.  Arc  (t,  s)  has  a  negatively  large  cost  and 
flow  below  its  upper  bound,  and  is  therefore  out  of  kilter,  while  ail 
other  arcs  are  in  kilter.  By  the  labeling  -’lies,  all  arcs  other  than 
(t,  s )  will  stay  in  kilter. 

We  first  label  node  s,  [t+,  u  -  x  j,  and  search  for  nodes  which 

ts  ts 

can  be  labeled  from  s.  There  are  initially  several,  since  c  .  =  0  and 

s  ■} 

xg  <  Ug  ^  on  all  arcs  (s,  j).  We  continue  the  labeling  process  until 
th-ire  is  a  breakthrough.  We  then  augment  flow  on  the  breakthrough  path 
b>  the  OKA  rules. 

Since  we  assigned  the  capacity  on  the  return  arc  (t,  s)  a  very 

VHr 

large  value,  the  arc  will  still  be  out  of  kilter.  Therefore,  repeat 

the  labeling  process  arid  further  breakthroughs  will  continue  to  augment 
flow  from  s  to  t.  Finally,  there  will  be  no  path  from  s  to  t  on  which 
every  arc  has  flow  less  than  capacity.  At  this  point  the  algorithm  will 

'k'k-k 

consider  all  arcs  connecting  labeled  to  unlabeled  nodes,  select  the 
minimum  cost,  c  ,  on  these  arcs,  add  this  value  to  all  node  prices, 
on  the  unlabeled  nodes,  and  attempt  to  laoel  them.  Generally,  a  label 
would  be  assigned  to  at  least  the  forward  node  of  the  arc  on  which  the 
minimum  c  had  been  found,  for  the  new  c.  or,  that  arc  would  be  zero. 

1  J  1  < 

In  this  case,  however,  flow  equals  capacity  on  all  arcs  connecting 
labeled  and  ur1 ah"led  nodes.  Consequently,  even  with  the  repricing, 
no  new  labels  will  he  assigned.  The  algorithm  will  once  again  consider 
all  arcs  connecting  labeled  to  nniabeled  nodes,  but  there  will  be  at 


‘if  a  solution  is  possible  with  positive  flow  (u,  .  0  on  some  path) 

:here  will  be  a  breakthrough ,  since  initially  t low  on  dll  arcs  is  zero. 


-It*  _ 

c 


0 


cs 


is  negativelv  very  large,  and  x 
!  S  r  S 


t  s 


These  arcs  will  now  have  x. 


i  1 
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least  one  less  arc,  since  bv  o^r  repricing,  at  least  one  c  now  equals 

zpro.  The  OKA  will  thus  continue  to  look  at  arcs  connecting  labeled  to 

unlabeled  nodes  which  have  c  ,  >  0,  it  will  increase  "  ,  and  continue  to 

ij  1 

remove  arcs  from  further  cons ide to 1 ion .  Eventually  the  only  arc  left 

★ 

for  consideration  with  a  labeled  and  unlabeled  node  will  be  (t,  s). 

At  this  point,  "  is  increased  by  a  sufficiently  large  number  so  that 
c  ^  =  0,  Arc  (s,  t )  is  now  in  kilter.  We  therefore  have  maximum  flow 
through  the  network. 

Notice  that  once  flow  augmentation  stops,  we  reach  a  maximum  flow. 
The  remainder  of  the  p.  cess  simply  computes  the  proper  prices.  We 
could  have  terminated  the  algorithm  as  soon  as  there  was  no  longer  a 
flow  augmenting  path  from  s  to  t. 

The  problem  of  finding  t  tie  shortest  route  through  an  uncapacitated 
network  is  identical  to  that  of  finding  the  minimum  cost  path  from  s  to 

t  through  a  network  on  which  each  arc  has  u  *  1,  and  on  which  a  return 

*■  J 

arc  (t  ,  s)  tias  t  =  0,  u  =1,  and  c  =  -M,  where  Mi?  v»*rv  large. 

’  '  ts  ts  ts  ’ 

The  OKA  will  send  one  unit  of  flow  from  s  to  t  at  minimum  cost.  Assume 

c.,  >  0  for  each  arc.  Start  the  problem  with  x,,  •  0  on  all  arcs,  and 
lj  ij 

”,  *  C  at  all  nodes.  Since  c  ,  >  0  on  all  arcs  except  arc  :t,  s),  the 
j  1J 

only  out -o  i-ki  It  e.  *■  arc  is  (t,  s).  By  the  labeling  rules,  node  s  is 


labeled  t  +  ,  e(s)',  where  e(s> 


ts 


ts 


fhe  algorithm  then 


attempts  to  label  nodes  from  s,  but  since  c,  >  0,  it  cannot.  It  there- 

l  j 

fore  revises  all  ”,  at  all  unlabeled  nodes  bv  adding  to  thorn  the  minimum 
J  , 

c ,  .  Cost  on  at  least  one  arc  (s  ,  J)  then  becomes  c  ■  0,  and  on 

i  J  5  ) 

this  arc  flow  is  less  than  capacity.  The  lorwnrd  node,  j,  of  this  arc 


receives  the  label  ^s+,  e(p',  where  c(j)  ■  min  e(s),  u 


Since 


ii  is  1,  e(j)  *  c  ( s )  ,  which  was  i.  This  indicates  that  the  second  label 
s  J 

at  a  1 1  labeled  nodes  will  be  one,  and  t he  problem  is  only  one  of  creaking 
through  on  some  pact,.  The  algorithm  continues  forward  from  all  labeled 
nodes  ,  finding  ail  c  '•  0,  and  being  forced  to  increase  ”  by  the  ap- 

1  .1  j 

propria  ft*  minimi, 


v  j 


This  then  alLovs  labeling  one  more  arc  to. ward. 


Eventually,  node  t  is  reached,  and  flow  along  the  breakthrough  path  is 


It  will  be  the  last  one  because  e 


verv  large.  Hence, 
** 


s  t 


M 


“l8. 


-M’  U5t 


Mi  ,  where  M  is 


This  is  the  feasible  but  nonoptimal  condition  described  on  page  la. 
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increased  by  one  unit.  All  arcs  along  the  path  are  in  kilter  since 


c  ■  0  for  each  arc  on  the  path,  and  x 
route  (or  minimum  cost  path)  we  need  only 


•*  u .  .  To 
ij 


trace  back 


find  the  shortest 
through  the  network 


from  t  to  s,  using  the  first  component  of  the  labels.  The  length  or 


cost  of  this  path  is 


given  by 
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